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ABSTRACT 


The  problem  we  address  is  that  of  choosing  a  deployment 
and  firing  doctrine  for  defending  separated  point  targets 
of  (potentially)  different  values  against  an  attack  by  an 
unknown  number  of  sequentially  arriving  missiles.  We  minimize 
the  total  number  of  defenders  subject  to  an  upper  bound  on  the 
maximum  expected  value  damage  per  attacking  weapon.  We  show 
that  the  Greedy  Algorithm  produces  an  optimal  integral  solu¬ 
tion  to  this  problem. 


v 


INTRODUCTION 


1  . 

In  the  late  1950s  and  early  1960s,  R.C.  Prim  and  W.T.  Read 
[5]  at  the  Bell  Telephone  Laboratories  and  the  Office  of  the 
Secretary  of  Defense  developed  a  method  to  deploy  and  fire  inter¬ 
ceptor  missiles  defending  a  collection  of  separated  point  targets 
against  an  attack  by  an  unknown  number  of  sequentially  arriving 
ballistic  missiles.  Basically,  the  scheme  involved  deploying 
and  firing  the  interceptors  in  a  manner  that  equalizes  the  prob¬ 
abilities  that  each  of  a  prescribed  number  of  attacking  weapons 
destroys  the  target.  These  deployments  (and  associated  firing 
doctrines)  have  become  known  as  "Prim-Read"  deployments. 

In  January  1964  a  major  study  of  limiting  damage  by  civil 
defense  and  antiballistic  missile  defense  was  completed  [1] . 

This  study,  performed  under  the  direction  of  Glenn  A.  Kent,  was 
highly  influential  in  shaping  United  States  policy  with  respect 
to  strategic  defense  (see,  for  instance,  [3]  for  background 
discussion) .  The  study  used  Prim-Read  deployments  and  firing 
doctrines  for  interceptors. 

Optimality  properties  of  the  Prim-Read  concept  were  estab¬ 
lished  by  Karr  in  [4] .  Some  background  discussion  of  the 
development  of  the  methodology  is  also  contained  in  that  docu¬ 
ment.  The  work  to  date  has  been  limited  to  continuous  (i.e., 
not  necessarily  integral)  firing  doctrines. 

In  the  present  paper,  we  focus  on  a  particular  class  of 
deployment  problems  whose  continuous  solutions  are  shown  in 
[4]  to  be  of  the  Prim-Read  variety,  and  develop  a  method  that 

produce  globally  optimal  solutions  of  integer  versions 
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of  these  problems.  This  class  will  thus  be  distinguished  by 
two  features:  (a)  the  number  of  defending  interceptors  allo¬ 
cated  to  each  attacking  missile  is  required  to  be  an  integer, 
and  (b)  these  problems  have  non-integer  versions  for  which 
(non-integer)  Prim-Read  deployments  are  known  to  be  solutions. 

The  underlying  "physics”  of  the  model,  as  described  in  the 
next  section,  is  that  of  [1,4]  and  [2,  p.  92].  The  optimization 
problem  that  we  study  is  taken  from  [4] ,  with  the  sole  modifi¬ 
cation  that  we  add  integer  restrictions  on  the  defense.  The 
authors  do  not  wish  to  imply  that  the  assumptions  here  are  the 
only  "correct"  ones;  numerous  alternative  models  can  reasonably 
be  addressed. 

The  problem  that  we  initially  address  has  the  general  form 
of  minimizing  the  total  number  of  interceptors  required  to 
defend  T  targets  against  an  attack  of  k  missiles,  where  neither 
k  nor  upper  or  lower  bounds  on  k  are  known  to  the  defender  in 

advance ,  subject  to  a  given  upper  bound  on  the  maximum  target 
value  destroyed  per  attacking  weapon.  In  the  next  section,  we 
will  introduce  the  specific  problem  assumptions  and  notation 
used  to  define  the  problem. 

In  Section  3,  we  note  the  existence  of  a  solution,  and 
show  that  the  multi-target  problem  reduces  to  a  collection  of 
single-target  problems.  Section  4  presents  some  properties 
of  solutions  to  the  single  target  problem. 

Section  5  contains  the  main  result  of  this  paper,  namely 
that  the  "Greedy  Algorithm"  solves  the  single-target  case. 
Combined  with  the  reduction  result  of  Section  3,  this  completely 
solves  the  specific  "Prim-Read"  problem  that  we  are  addressing. 

Section  6  treats  the  special  case  of  perfect  interceptors, 
whose  solution  is  very  easily  obtained  and  is  shown  also  to  be 
a  limiting  case  of  the  solutions  produced  by  the  algorithm  of 
Section  5.  Section  7  presents  two  detailed  examples. 
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ASSUMPTIONS  AND  NOTATION 


2  . 

The  general  problem  we  address  has  the  form  of  minimizing 
the  total  number  of  interceptors  required  to  ensure  that  the 
maximum  possible  damage  under  any  attack  is  bounded  by  a  linear 
function  of  the  attack  size.  The  purpose  of  this  section  is  to 
introduce  specific  assumptions  and  notation  so  that  the  problem 
becomes  well-defined. 

(1)  We  assume  that  the  defender  must  set  his  strategy 
first,  and  in  ignorance  of  both  the  attacker's  resources  to  be 
expended  (i.e.,  total  number  of  attacking  weapons)  and  inten¬ 
tions  concerning  allocation  of  missiles  among  targets  (i.e.,  the 
attacker's  strategy).  We  allow  the  attacker  knowledge  of  the 
specific  defense  strategy  adopted  by  the  defender. 

(2)  We  assume  that  there  are  T  targets,  and  that  each 
target  i  has  a  value  v(i)  >  0.  The  v(i)  need  not  be  integers. 

(3)  We  assume  that  defending  interceptors  must  be  assigned 
to  specific  targets,  i.e.,  they  cannot  be  used  in  an  "area 
defense"  mode.  We  also  assume  that  the  targets  are  separated 
to  the  extent  that  an  interceptor  assigned  to  defend  a  speci¬ 
fic  target  cannot  be  used  to  defend  a  different  target.  More¬ 
over,  an  attacking  missile  directed  at  a  specific  target  has 

no  effect  on  the  survival  or  destruction  of  a  different 
target . 

(4)  We  assume  that  attacking  weapons  arrive  sequentially 
(i.e.,  so  that  they  can  be  labeled  "first",  "second",  ...)  with 
sufficient  time  between  arrivals  that  the  fate  of  any  particu¬ 
lar  attacking  missile  is  independent  of  that  of  every  other 
attacking  missile.  The  defender  does  not  know  either  the 
attack  size  or  any  (nontrivial)  upper  or  lower  bounds  on  it. 

The  number  of  interceptors  committed  to  a  given  attacking 
missile  must  be  chosen  with  no  knowledge  concerning  how  many 
missiles  will  follow  at  the  same  target.  In  effect  the 


defender  is  thereby  forced  to  assign  interceptors  in  advance  to 
potentially  arriving  missiles  at  each  target.  • 

(5)  We  assume  that  neither  side  can  change  his  strategy 
during  the  attack.  Thus  we  are  considering  only  non-adaptive 
strategies  for  both  sides. 

Based  on  the  above  assumptions,  we  define  a  defense 
strategy  d  to  be  a  semi-infinite  matrix 


d  = 


d( 1 , 1) 

d( 1, 2)  ... 

d(i,l) 

d( i , 2 )  ... 

d(T,l) 

d(T, 2)  ... 

where  d(i,j)  -  number  of  interceptors  assigned  at  target  i,  to 

t/ h 

be  directed  at  the  j  incoming  attacking  weapon  (if  there  is 

one)  at  that  target.  We  will  also  use  the  notation  d(i)  to 
th 

denote  the  i  row  of  d.  Throughout  this  paper,  the  d(i,j)'s 
are  restricted  to  be  nonnegative  integers .  In  practice,  all 
but  a  finite  number  of  the  d(i,j)'s  will  be  zero. 

In  a  similar  fashion,  an  attack  strategy  is  a  vector 


a  = 


a(  1) 


a(i) 


a(T) 


4 


where  a(i)  denotes  the  number  of  attacking  weapons  to  be 
directed  at  target  i;  the  a(i)'s  are  also  restricted  to  be 
nonnegative  integers. 

We  remind  the  reader  that  the  defender  must  choose  d 
first,  and  the  attacker  is  then  free  to  choose  a,  based  on 
knowledge  of  d.  Note  that  we  are  giving  the  attacker  knowledge 
of  each  d(i,j),  i.e.,  he  not  only  knows  the  number  of  inter¬ 
ceptors  stationed  at  target  i,  but  he  also  knows  the  schedule 
by  which  these  interceptors  will  be  fired. 

(6)  Let  q  t  [0,1)  denote  the  probability  that  an  inter¬ 
ceptor  engaging  an  attacking  missile  fails  to  destroy  that 
missile.  We  assume  that  different  interceptors  deployed 
against  a  given  attacking  missile  are  independent,  and  that  the 
same  q  is  valid  for  every  interceptor  and  attacking  missile 
engagement  no  matter  how  many  interceptors  attack  each  missile. 

(7)  We  assume  that  an  attacking  missile  that  penetrates 
the  defense  will  destroy  its  intended  target  with  probability 
one . 


For  a  given  deployment  d  and  a  given  attack  strategy  a, 
the  probability  that  target  i  is  destroyed  is 

a(i ) 

n 

j=l 


1  -  n  (l-qd(l’J'  )  )  • 


0 


with  the  convention  that  n  (•)  =1,  this  is  zero  if  a(i)  =  0. 


Therefore 


1 


T 


V(d,a)  =  l  v(i) 
1=1 


a  ( i )  .  / .  .  v 

1  -  It  ( l-qd  1  ^  ) 

j=l 


It  is  this  quantity  that  an  attacker  might  seek  to  maximize  for 
a  given  d,  subject  to  an  upper  bound  on  the  number  of  missiles 
to  be  expended. 


If  the  number  of  attacker's  missiles  to  be  expended  were 
known  to  the  defender  in  advance  to  be,  say,  k,  a  reasonable 
problem  to  address  might  be 

co 

minimize  [  J  d(i,j) 
d( i, j )  i=l  j  =  l 

subject  to  V(d,a)  <_  u, 

for  all  a  such  that 
T 

l  a(i)  <_  k  , 

i=l 

where  u  is  some  prescribed  upper  limit  on  the  expected  target 
value  damage.  However,  in  the  current  context,  k  is  not  known 
In  advance.  The  alternative  that  we  address  assumes  instead  a 
bound  on  the  maximum  expected  damage  incurred  per  attacking 
weapon .  This  alternative  is  in  the  same  spirit,  but  makes  sense 
when  k  is  not  known  in  advance. 

3.  THE  MULTI-TARGET  PROBLEM  AND  SOME  OF  ITS  PROPERTIES 

In  Sections  3-5  we  assume  that  0  •<  q  <  1 .  The  case  q  =  0, 
in  which  a  single  interceptor  is  certain  to  destroy  an  attack¬ 
ing  missile,  is_  of  physical  interest  and  is  treated  in  Section 
6  using  arguments  different  from  those  in  Sections  3-5.  As  will 
emerge,  the  case  q = 0  is  much  easier  to  handle  because  non¬ 
integral  numbers  of  interceptors  would  never  arise  in  optimal 
solutions  even  if  they  were  not  forbidden,  yet  still  this  case 
can  be  regarded  as  a  limiting  verion  of  the  case  q>0. 

Let  a  real  number  s  >  0  be  given.  This  number  will  serve 
as  an  upper  bound  on  the  maximum  expected  damage  incurred  per 
attacking  weapon.  Some  schemes  for  selecting  reasonable  values 
of  s  are  given  in  [3];  we  note  just  one  here. 

If  V  is  the  sum  of  all  of  the  target  values,  then  V  is 
a  given  upper  bound  on  the  total  expected  damage  incurred,  and 
If  the  total  number  of  attacker's  weapons  equals  K,  then 
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s=V/K  is  the  corresponding  upper  bound  on  the  maximum  expected 
damage  incurred  per  attacking  weapon  is  the  attacker  commits 
his  entire  inventory.  The  problems  presented  below  ensure  that 
this  average  is  not  exceeded  even  if  the  attacker  holds  back 
some  missiles  or  builds  new  ones. 


The  problem  we  address  herein  has  the  form 

oo 

minimize  E  Ed(i,j) 
d(i,j)  i=l  j=l 

T 

subject  to  V( d, a)  £  s  XI  a(i),  for  all  a, 

i=l 

where ,  we  recall,  each  d(i,j)  and  each  a(i)  is  further 
restricted  to  be  a  non-negative  integer. 

We  note  that  problem  P(s)  has  the  equivalent  formulation 

00 

minimize  E  Ed(i,j) 
d ( i , j )  i=l  j=l 

T 

V(d,a):  E  a(i)<k  <  sk 

i=l  ) 

for  k  =  1,2,3,  ...  , 


subject  to  max 


P(s) 


and  we  will  address  this  form  in  several  of  the  proofs. 

#  Theorem  1  (Existence  of  a  Solution):  Problem  P(s)  has  a  solu¬ 

tion  for  each  value  of  s. 

Proof :  First  note  that  the  number  of  potentially  effective 


constraints  is  finite.  For1  if  we  set  k  = 


‘  T 

E  v(i)/s 

i=l 


we 


x0n  those  occasions  where  it  is  necessary  to  round  non-integer  quantities, 
we  use  the  notation 


llxjj  =  max{i:i  is  integer,  i  <  x} 
fix])  =  min{i:i  is  integer,  i  >  x}. 
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note  that  the  constraints  for  k,  k  +  1,  ...  are  irrelevant, 

T  • 

since  V(d,a)  is  bounded  above  by  X]  v(i). 

i=l 


Let 


f,(d)  =  max'jvC  d,  a)  :  XT  a(  i)  <kl 

1  i=l  ’ 


denote  the  kth  constraining  function,  so  that  the  kth  constraint 

has  the  form  • 


fk(d)  <  sk. 

Now  define  the  function  • 


gk(d)  =  v  l  (1-  n  (l-qd(ljJ))) 

i=l  j=l 

where  v  =  max{v( i) : i-1, . . . ,T} .  Since  fk(d)  £  gk(d),  it  suffices 
to  establish  feasibility  for  the  constraints  gk(d)  <  sk 

(k=l, . . . ,k-l) .  Prom  here,  it  is  easy  to  show  the  existence 
of  feasible  solutions  of  the  form 


D  for  j  =  1,...,  k-1,  i  =  1,...,  T 
0  for  j  =  k,  .  .  .  ;  i  =  1,...,  T  , 


where  D  is  chosen  large  enough  that  <  ( 1—  ( 1  -  —  )~^) 

Tv 

for  k=l, . . . ,  k=l,  so  that  P(s)  is  feasible.  This,  together  with 
the  constraints  that  d(i,j)  take  on  integer  values,  and  the  form 
of  the  objective  function,  guarantees  that  P(s)  has  a  solution.  D 

The  next  theorem  establishes  the  result  that  the  problem 
PCs)  can  be  solved  by  solving  a  collection  of  single-target 
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defense  problems,  one  for  each  target.  This  reduction  permits 
us  to  focus  attention  in  Sections  4  and  5  on  the  single-target 
problem  P(r)  to  be  stated  momentarily. 

For  convenience  of  notation,  we  define 


V. (d(i) 

l 


k)  =  v(i)  (1  - 


(> 


k 

n  ( l-q 
0=1 


d  ( i ,  j  ) 


)■ 


where  d(i)  =  (d(i,l) ,d(i,2) , . . . )  as  above.  Then 


T 

V(d,a)  =  l  V. (d(i) ,a(i) )  . 
i=l  1 


The  proof  of  the  following  theorem  depends  on  the  fact  that  the 
expected  damage  is  the  sum  of  the  expected  damage  at  each  target, 
and  this  is  bounded  above  by  the  function  sk,  which  is  linear 
in  k . 

Theorem  2  (Reduction  of  the  Multi-target  Problem) :  For  each 
i  e  {1,...,T},  let  d*(i)  solve  the  problem 

oo 

(minimize  £  d(i,j) 

0=1 

subject  to  V.(d(i),k)  j<  sk  (k=l,2,...). 

Then  d*  =  (d* ( 1 ), d* ( 2 ),..., d* (T ) )  solves  P(s). 

Proof:  First  note  that  d*  is  feasible  for  problem  P(s),  since 

V.  (d*  (i)  ,a(I ) )  <_  sad) 

for  any  a(i),  so  that 

maxjv(d*,a): J^aCi)  <  k|  <  sk 
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for  each  k  =  1,  2 ,  ...  . 

Now  suppose  that  d*  is  not  optimal  for  P(s);  then  by 
Theorem  1  there  exists  a  feasible  deployment  d  which  satisfies 

T  oo  t  00 

l  l  d(i,j )  <  l  l  d* ( i j j )  , 

i=l  j=l  i=l  j=l 

so  that,  at  some  target  I,  we  have 


I 


j  =  l 


d(l,j) 


<  I  d»(i,j)  . 

j=l 


By  definition. 


V  (d(I),k)  =  l  V . ( d( i ) , 0 )  +  VT(d(I),k), 
i=l  1  1 

ItO: 

=  V(d,k) 


where  k  is  the  Ith  unit  vector  times  k.  Since 


V(d,k)  <  maxjv(d,a) :  J  a(i)<kj 


<  sk  , 


we  see  that  d(I)  is  feasible  to  problem  P_(s),  but  entails  a 
smaller  total  deployment  than  d*(I),  contradicting  the  assumed 
optimality  of  the  latter  for  P  (s).  [] 

Since  we  are  now  able  to  focus  attention  on  the  single¬ 
target  problem,  it  is  convenient  to  simplify  the  notation.  To 
that  end,  we  will  now  employ  the  symbol  d  to  denote  a  vector 

whose  j _ component  d(,j)  denotes  the  number  of  interceptors 

assigned  to  attacking  weapon  j.  Also,  we  introduce  the  symbol 
r  in  place  of  s/v(i)  and  suppress  reference  to  target  i.  Note 
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that  the  values  of  interest  of  r  lie  between  zero  and  one  (the 
value  of  a  single  target  is  normalized  to  equal  one) . 


The  reader  should  note  that  this  change  of  notation 
(previously  d  was  a  matrix ,  and  d(i)  was  a  vector)  will  remain 
in  place  throughout  the  rest  of  the  paper,  except  in  Section 
6  and  Example  2  of  Section  7,  where  a  multi-target  problem  is 
solved . 


The  single-target  problem  thus  becomes 


minimize  H  d(j) 
P  (  r )  :  J'_1 


subject  to  V(d,k)  £  rk  (k=l,2,...)  , 


where  V(d,k) 


£  (l-qaU>)J 


4.  THE  SINGLE-TARGET  CAS E- -M0N0T0 N I C I TY  AND  ADMISSIBILITY 

Definition:  Let  d  denote  any  deployment  vector  ( d( 1) , d( 2 ),...) . 

We  say  that  d  is  monotone  if 

d(  1)  >  d(  2)  >  .  .  .  . 


Under  an  attack  by  missiles  arriving  sequentially  in  time, 
a  monotone  deployment  defends  more  heavily  against  the  earlier 
attacking  weapons,  an  action  which  seems  reasonable. 

During  the  proof  of  Theorem  3  and  for  the  remainder  of  the 
paper  we  suppress  dependence  on  r  and  write  "Problem  P(r)"  as 
simply  "Problem  P" . 

Theorem  3:  For  every  r  there  exists  a  monotone  optimal  solution 
of  problem  P. 
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Proof:  Let  d*  denote  an  optimal  solution  of  problem  P. 


Then 


V(d*,k)  =  1  -  n  (l-qd*(j'))  <  rkj  k=l, 2 , . . .  . 

j=l 

If  d*(j)  <  d*(j+l)  for  some  j ,  define  the  vector  d  from  d*  by 
interchanging  the  components  d*(J)  and  d*(J+l).  Both  d  and 
d*  have  the  same  component  sum,  and  we  note  that 

V(d,k)  =  V(d*,k)  for  k  ^  J 

while 

V(d,J)  <  V(d*J) 

so  that  d  is  also  feasible  to  problem  P.  Since  it  is  easily 
seen  that  one  can  repeatedly  interchange  adjacent  components 
of  d*  that  increase,  until  one  obtains  a  monotone  solution  of 
P,  the  theorem  is  proven.  0 

Note  that  the  above  interchange  strictly  decreased  the 
value  of  V(d*,J),  while  keeping  the  other  quantities  V(d*,k) 
unchanged,  i.e.,  with  the  same  total  number  of  interceptors, 
the  defender  can  limit  his  expected  damage  under  an  attack  by 
k  ^  J  missiles  to  the  same  amount,  and  decrease  this  damage  if 
attacked  by  J  missiles.  This  is  obviously  a  desirable  property 
from  the  defender's  point  of  view,  and  motivates  us  to  make  the 
following  definition. 

D e f i n  i  t  i  o  n  :  A  deployment  d  that  solves  problem  P  is  admissible 
if  there  is  no  other  deployment  d  such  that 

(a)  d  also  solves  problem  P 

(b)  V(d,k)  <_  V(d,k)  for  all  k 

and  V(d,k) <  V(d,k)  for  some  k  . 


The  proof  of  the  previous  theorem  establishes  that  admis¬ 
sible  deployments  are  monotone  (or  else  they  could  be  "monoton- 
ized",  yielding  deployments  d  having  at  least  one  strictly 
better  value  V(d,k)).  The  next  theorem  establishes  that 
optimal  monotone  deployments  are,  in  fact,  admissible. 

Theorem  4 :  An  optimal  deployment  dR  of  problem  P  is  monotone 
if  and  only  if  it  is  admissible. 

Proof:  Let  d*  be  an  optimal  monotone  deployment,  and  assume 

that  it  is  not  admissible.  Thus  there  must  be  another  optimal 
deployment  d  that  also  solves  P,  and  for  which 

V  (  d ,  k )  <  V(d*,k) 

for  all  k,  with  strict  inequality  holding  for  some  k.  We  can 
assume  that  d  is  also  monotone  (perhaps  "monotonized"  by  the 
procedure  established  above). 

Since  V(d,l)  £  V(d*,l),  it  follows  that  d(l)  >_  d*(l). 

Let  k-^  denote  the  smallest  value  of  k  such  that  d(k^)  >  d*(k-^). 
We  note,  then,  that 


V(d,k1)  <  V(.d*,k1)  . 

00 

Because  both  d  and  d*  are  optimal,  we  know  that  E  dtk)  = 

oo  k=l 

E  d*(k).  Now  let  k0  denote  the  smallest  value  of  k  >  k.  such 
k=l  d  1 

that 


2^  2 
l  d (k )  =  l  d*(k). 
k=l  k=l 


We  claim  that  V(d,k2)  >V(d*,k2) 
that  d*  was  not  admissible.  It 


which  violates  the  assumption 
suffices  to  show  that 


n2(iV(k))  <  n2  (i-qd*(k)) 

k=l  k=l 


and  this  is  a  direct  result  of  the  following  Lemma,  if  we  set 
t  =  k2  -  k1  +  1,  ai  =  dCi+k-,-1),  and  b±  =  d«(i+k1-l).  0 

Lemma :  Let  a  =  (a^)  and  b  =  ( )  be  a  pair  of  positive 

monotonically  nonincreasing  vectors  of  size  t  such  that 


(a)  E  a,  >  E  b . 

i=l  1  i=l  1 


for  s  =  1,2, ,  t-1  ; 


(b )  E  a,  =  E  ( 1-q  i)  . 


i=l 


i=l 


Then  for  q  e (0 ,1 ) , 


t  a.  t  b. 

II  (1-q  x)  <  II  (1-q  x) 
i=l  i=l 


Proof:  See  Appendix  A. 


There  may  be  more  than  one  optimal  admissible  solution,  as  • 

the  following  example  exhibits. 

< 

Example:  With  r  =  1/3  and  q  =  0;6,  the  deployments  d*  = 

(3,2)  and  3  =  (4,1)  are  both  optimal,  admissible  solutions.  ^ 

We  have 


V(d*,l)  =  .2160 
V(d*,2)  =  .4982 
V(d*,3)  =  1.000 


V( d, 1 )  =  .1296 
V(d, 2)  =  .6518 
V(d,3)  =  1.000  , 
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so  that  the  solution  (3,2)  yields  a  higher  expected  damage  if 
there  is  but  a  single  attacker,  but  a  lower  expected  damage  if 
there  are  two  attackers .  The  average  expected  damage  per 
attacker  values  are 


V  (  d  *  ,  1 ) 
1 

V(d* , 2) 
2 

V(d*,3) 


3 


0.2160 

0.2491 

0.3333 


V  (.  d ,  1 ) 
1 

V(d,2) 

2 

V(  d , 3) 

3 


0.1296 
0.3259 
0.3333  , 


which  are  all,  as  required,  bounded  above  by  r  =  1/3- 


5.  SOLUTION  OF  THE  S I NGL E - TAR  GET  PROBLEM 

In  this  section,  we  present  a  proof  that  the  "Greedy 
Algorithm"  solves  the  single-target  problem.  The  method  begins 
by  computing  the  minimum  number  of  defending  interceptors 
necessary  to  ensure  that  the  expected  value  of  the  damage  from 
one  attacker  is  at  most  r.  The  iteration  step  assumes  that 
k-1  attacking  weapons  have  been  assigned  interceptors,  and 
the  minimum  number  of  interceptors  to  be  assigned  to  attacking 
weapon  k  is  then  computed,  subject  to  the  usual  bound. 


The  Algorithm 


Let  cT(l)  = 


Jin  r 


Jin  q 


and,  given  d(l) , . . . , 


d(k)  = 


(1- 


1-rk 


*nl  kiMl-qd(l)) 

i=l _ 

In  q 


d(k-l) ,  let 
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for  2  <_  k  <  1/r.  For  k  _>  1/r,  let  d(k)  =  0.  (Note  that  the 

general  form  for  k  >  2  is  actually  valid  for  k  =  1,  with  the 

0 

convention  that  n  (•)  =  1.) 

i=l 

Theorem  5 :  The  above  recursion  yields  an  optimal  solution  of 
d  of  Problem  P. 


Proof:  The  first  component  d(l)  solves 


and  the  k^d 


minimize  d(l) 

subject  to  1  -  (l-qd  d  )  £  r 
component  d(k)  solves  the  problem 

minimize  d(k) 
subject  to 


/k-1 

-(  n  1 

\£=1 


■) 


(l-qdU)))  (l-qd(k))  <  rk 


f°^  2  <k  <  1/r,  as  is  easily  verified. 

Let  n  =  |J_l/rjJ  .  Clearly,  if  d(l  ),  ...,  d(n)  is  any  sequence 
feasible  for  the  first  n  constraints,  then  d(l),  ...,  d(n),0,0, 
...  is  feasible  for  all  constraints,  as  the  (n+l)s1:,  (n+2)nd, 

...  are  redundant. 

Let 


C  =  { x ( Rn : x£0 ,  V ( x , k )  £  rk ;  k= 1 , . . . , n } 

(the  members  of  C  are  not  restricted  to  have  integer  components), 
and  let  L  denote  the  lattice  of  points  with  integer  coordinates 
in  (Rn)+.  We  note  that 


d  «  C  0  L 


3 
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and  we  must  show  that  d  e  COL  implies  that 


n 


n 


l  d(i)  >  l  d(i)  . 


i=l 


i=l 


We  define,  for  each  £  =  1 


o 

.  n,  a  deployment  d  by 


setting 


d  ( i ) 


if  i  /  i 


d( £)  -  1  if  i  =  £ 


None  of  these  deployments  can  be  feasible,  because  d(£)  was 
the  smallest  d(£)  to  satisfy  constraint  £,  given  d(l), 
d(£-l) . 


Let 


n  n 


K  =  (d:  l  d(i)  =  l  d(i)  -  1} 
1=1  i=l 


£ 

so  that  we  need  to  show  that  CO  (KHL)  =  0.  The  points  d 


are  those  elements  in  KOL  that  are  "nearest"  to  d.  In  what 

follows,  we  will  show  that  other  points  in  KHL  (which  are 

farther  away  from  d)  would  have  more  difficulty  in  being 

£ 

feasible  than  any  of  the  d  's,  so  that  all  of  KHL  lies  outside 
of  C. 

Any  element  d  in  KOL  can  be  obtained  from  an  element  in 

{d"L,...,dn}  by  a  sequence  of  transformations  that  simultaneously 

£ 

increase  a  component  of  d  by  1,  and  decrease  another  component 
of  d£  by  1. 

Now  fix  some  £,  and  let  i,  and  ip  be  distinct  integers 
between  1  and  n.  Define 


|d£(i1)  +1  i  =  i1 

dA(i2)  -  1  i  =  ±2 

d£(i)  i  7"  i^ ,  i 2  , 

so  that  d  is  in  KfiL.  We  want  to  show  that  d  /  C;  for  other¬ 
wise,  d  would  be  a  better  feasible  solution  than  d. 

The  first  component  of  d  that  differs  from  d  must  offer 
an  increase  in  that  component  because  of  the  way  d  was  con¬ 
structed.  But  we  can  assume  that  d  differs  from  d  in  components 
i-p  and  i2  (if  i^  =  i,  then  d  =  d  ,  which  is  known  to  lie 
outside  of  C) .  Thus,  we  can  assume  that 

i^  <  min{i2,£}  . 


We  claim  that 


V(d,m)  _>  V( d£ ,m)  ,  (1) 

where  m  =  max  {i2,£}.  Since  we  know  that  V(d£,m)  >  rm,  it 
will  follow  that  d  /  C. 

Rewriting  (1),  we  see  that  it  is  equivalent  to 


d V.  )+l  d£(i  )-l  d£(i,)  d£(i„) 

( 1-q  )  ( l-q  )  <  (1-q  )  ( 1-q  2  )  ,. 


which  further  is  equivalent  to 


d£(in)  d£(  ip ) 

q  1  (1-q)  <  q  (q  -1)  , 
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i .  e .  , 


d£(i1 )  -  d£(i?)+l 

q  <  1  , 

which  is  certainly  true  if  d^d^)  -  d£(i2)  +  1  0  ,  i.e.,  if 

d( i1)  +  1  _>  d(  i2)  (2) 

(d£  is  the  same  as  d,  except  in  component  Z)  .  Now,  while  the 
d  produced  by  the  Greedy  Algorithm  may  not  be  monotone,  at 
least  we  see  by  construction  that  we  always  have  d(i)  +  1  > 
d(i+l),  so  that  (2)  is  satisfied. 

The  procedure  extends  inductively  to  any  element  in  KOL 
and  the  theorem  is  proved.  [j 

Before  proceeding  to  the  next  section,  we  point  out 
that  Burr  has  derived  a  different  algorithm  based  on  an 
equivalent  formulation  of  the  problem.  It  is  included  in 
Appendix  B  of  this  paper,  and  has  a  similar  "greedy"  nature, 
but  (unlike  this  method)  always  produces  monotone  deployments. 

6.  THE  CASE  q  =  0 

For  this  section  we  revert  to  the  multiple-target  formula¬ 
tion  and  notation  of  Sections  2  and  3.  When  q  =  0 ,  interceptors 
are  perfectly  effective  and  the  defender  ensures  destruction  of 
an  attacking  missile  by  committing  just  one  interceptor  against 
it.  Hence  we  may  and  do  restrict  attention  to  deployments 
d  =  (d(i,j))  for  which  d(i,j)  is  either  zero  or  one  for  each  i 
and  j.  Furthermore,  it  is  evident  that  no  optimal  solution  d* 
to  the  basic  problem  P(s)  can  have  (for  some  i)  d*(i,j)  =  0  but 
d*(i,jO  =  1  for  some  j  <  j'  (since  the  deployment  obtained  by 
changing  d*(i,j")  to  zero  has  the  same  payoff  function  but 
strictly  fewer  interceptors).  Thus  we  consider  only  deployments 
d  of  the  form 
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j=l, . . . ,  S(i) 

J  >  5  (  i )  , 

where  the  S(i)  are  nonnegative  integers.  (If  6 ( i )  =  0,  then 
d(i,j)  =  0  for  all  j.)  The  deployment  d  is  completely  speci¬ 
fied  by  the  6(i),  so  we  take  the  6(i)  as  the  choice  variables 
of  the  problem,  which  now  becomes 


1 } 

d  ( i , j  )  =| 

0, 


P'(s): 


T 

minimize  j  6(i) 

6(i)  1=1 

T 

subject  to  V(d,a)  <  s  J  a(i), 

i=l 


for  all  a  . 


The  interpretation  of  <5(i)  is  the  number  of  attacking 

missiles  against  which  target  i  is  defended  and,  because  q  =  0, 

s  t 

to  which  it  is  invulnerable.  Against  the  (5(i)+l)  attacking 
missile  it  is  completely  defenseless  and  would  be  destroyed.  It 
follows  that  the  payoff  function  V  assumes  the  simplified  form 

T 

V(d,a)  =  v(i)I(a(i)>6(i) )  , 

i=l 


where  for  an  event  A,  1(A)  is  the  indicator  function  of  A.  (If 
A  occurs,  then  1(A)  =  1  and  otherwise  1(A)  =  0.) 

Problem  P'(s)  is  virtually  trivial  to  solve. 

Theorem  6:  The  unique  solution  to  Problem  P"*(s)  is  given  by 

6*(i)  =  ,  1=1,..., T  . 

Proof:  If  the  constraints  of  P^(s)  are  satisfied,  then  in 

particular 


# 
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v(i)  £  s(<5(i)+l)  (3) 

for  every  i.  Conversely,  If  (3)  holds  for  every  1,  then  for 
each  attack  allocation  a 

T 

V(d,a)  <  E  s(6(i)+l)I(a(i)>6(i)+l) 

1=1 

T 

£  s  Ea(i)  . 

1=1 

Thus  Problem  P^(s)  reduces  to  the  problem 

T 

minimize  E  6(1) 

6(1)  1-1. 

subject  to  v ( 1 )  £  s(5(i)+l)  ,  1=1,..., T  , 

and  of  course,  subject  also  to  the  constraint  that  each  6(1) 
be  a  nonnegative  integer.  The  unique  optimal  solution  to  this 
latter  problem  is  as  given  above.  [] 

The  optimal  deployment  6*  is  monotone  at  each  target, 
which  is  consistent  with  results  in  Section  4.  Also,  taking 
limits  as  q  -*  0  in  the  algorithm  of  Section  5  yields  the  deploy¬ 
ment 


d(k) 


lim 

q->0 


log(rk) 

log  q 


II,  if  k<v/s 

0,  if  kvv/s  , 

*  ^ 

which  is  the  same  deployment  as  6  from  Theorem  6. 


7.  EXAMPLES 

In  this  section,  we  present  two  examples  in  order  to 
illustrate  the  method  and  to  compare  the  results  with  those 
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obtained  by  ignoring  the  integer  restrictions.  The  first  exam¬ 
ple  addresses  a  single  target  case  and  is  solved  for  a  wide 
range  of  r  values,  while  the  second  example  addresses  a  nation¬ 
wide  case  for  one  value  of  the  parameter  s,  and  the  target 
values  set  to  be  the  populations  of  U.S.  cities. 


Solution  of  a  Single-Target  Case 

The  defender  wishes  to  defend  a  single  target  against  an 
attack  of  unknown  size.  He  decides  that  he  will  surrender  the 
target  if  8  warheads  are  directed  against  it,  but  wishes  to 
structure  a  defense  against  attack  sizes  of  0,1,..., 7.  The 
defender  decides  to  set  his  defense  so  that  the  maximum  proba¬ 
bility  of  target  destruction  (as  a  function  of  the  attack  size) 
is  bounded  above  by  a  linear  function  of  the  attack  size  k. 

Thus  he  decides  to  set  an  upper  bound  of  min  {k/8,l}  on  the 
probability  of  target  destruction. 


In  this  example,  the  defender's  interceptors  have  a  kill 
probability  of  (1-q)  =  0.4. 


According  to  the  algorithm,  we  compute 


d(  1) 


d(  2) 


Jin  0.125 
Jin  0 . 6 


Jin  1 


_  1  ~  2/8\ 

l-(0 . 6)  -V 

j£n  0 .  £ 


4 


d(3)  =  3 
d( 4)  =  4 

d( 5 )  =  3 
d(6)  =  2 
d(7)  =  2, 
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and  d ( 8 )  =  d(9)  =  •••  =0.  With  this  deployment,  the  defender 
will  assign  five  interceptors  to  the  first  attacking  missile, 
four  to  the  second,  and  so  on. 

Note  that  this  deployment  is  not  monotone.  An  equivalent 
monotone  solution  is  d  =  (5, 4, 4, 3, 3, 2, 2,0, . . . ) .  Both  d  and  d 
use  23  interceptors. 

The  attacker's  optimal  probability  of  target  destruction 
as  a  function  of  the  attack  size  is  listed  in  Table  1,  for  both 
d  and  d,  and  plotted  in  Figure  1.  Note  that  d  yields  a  better 
defense  against  an  attack  of  size  three  than  does  d. 

We  go  on  to  compute  the  optimal  defense  for  all  values  of 
r  in  the  range  [1/10,  1).  Table  2  contains  the  results.  In 
[4 ] ,  Karr  derives  the  optimal  total 


when  1/r  is  an  integer  and  the  integer  restrictions  on  the 
components  of  d  are  dropped.  These  totals  are  included  in 
Table  3-  Figure  2  displays  information  from  Tables  2  and  3- 

Solution  of  a  Nationwide  Case 

In  this  example,  we  take  cities  in  the  U.S.  as  targets, 
with  their  populations  in  thousands  as  values.  Population  data 
are  from  the  1980  Census  as  reported  in  the  1983  U.S.  Almanac. 
The  value  of  s  is  200  (the  defense  is  to  be  set  so  that  no 
single  attacking  weapon  can  kill  more  than  200,000  people)  and 
q  =  0.5-  Table  4  presents  the  continuous  solution  computed 
from  the  method  of  [4] ,  the  continuous  solution  rounded  up  for 
each  attacking  weapon  at  each  target,  and  the  integer  solution 
from  applying  the  algorithm  of  Section  4.  The  required  inter¬ 
ceptor  stockpiles  are  349,  446  and  4l4,  respectively,  as  shown 
in  Table  4. 


Table  1.  PROBABILITY  OF  TARGET  DESTRUCTION 


k 

V(cT,k) 

v(3,k) 

1 

0.07776 

0.07776 

2 

0.1 9728 

0.1  9728 

3 

0.37067 

0.30121 

4 

0.45223 

0.45223 

5 

0.57055 

0.57055 

6 

0.7251 5 

0.7251 5 

7 

0.8241 0 

0.8241  0 

8 

1  .0 

1  .0 

9 

1  .0 

1  .0 

I 
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PROBABILITIES  OF  TARGET  DESTRUCTION  VS  ATTACK  SIZE 


Table  2.  SOLUTION  OF  THE  SINGLE  TARGET  CASE 


_  v  Optimal 

r  Range  d  via  the  Greedy  Algorithm  Sum 
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Table  3.  SOLUTION  OF  THE  SINGLE  TARGET  CASE  WITHOUT  INTEGER  RESTRICTIONS 


27 


J  dO  NOLLONfld  V  SV  lN3WA01d30  1V101  =  (J)0 


CVJ 


a> 

s- 

=3 

CD 
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TOTAL  OPTIMAL  DEPLOYMENT 


Table  4.  SOLUTION  OF  A  NATIONWIDE  CASE 


Target 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 


Prim- 

Read 

Integer 

Stockpi 1 e 

Rounded 

Stockpil e 

134.92 

153 

139 

40.37 

49 

43 

40.20 

45 

42 

16.85 

21 

19 

15.22 

17 

16 

9.54 

14 

12 

5.93 

8 

8 

5.60 

8 

7 

4.48 

5 

5 

4.47 

5 

5 

4.30 

5 

5 

4.28 

5 

5 

3.72 

5 

5 

3.50 

5 

5 

3.15 

5 

5 

3.06 

5 

4 

3.04 

5 

4 

3.03 

5 

4 

2.42 

3 

3 

2.37 

3 

3 

2.35 

3 

3 

2.31 

3 

3 

2.21 

3 

3 

1  .86 

3 

3 

1.84 

3 

3 

1.55 

3 

3 

1  .52 

3 

3 

1  .47 

3 

3 

1.26 

3 

3 

1.26 

3 

3 

1.25 

3 

3 

1 .03 

3 

3 

0.94 

1 

1 

0.94 

1 

1 

0.89 

1 

1 

0.87 

1 

1 

0.85 

1 

1 

0.85 

1 

1 

0.84 

1 

1 

0.83 

1 

1 

0.79 

1 

1 

0.79 

1 

1 

0.76 

1 

1 

0.73 

1 

1 

0.73 

1 

1 

Val  ue 

7071 

3005 

2997 

1688 

1594 

1203 

904 

876 

787 

785 

765 

763 

701 

679 

646 

638 

637 

636 

574 

565 

563 

557 

541 

494 

491 

456 

453 

448 

425 

425 

424 

403 

385 

385 

371 

366 

361 

361 

358 

355 

347 

345 

339 

332 

331 
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Table  4. 


SOLUTION  OF  A  NATIONWIDE  CASE  (Continued) 


Target 

Val  ue 

Prim- 

Stockpile 

Read 

Rounded 

Integer 

Stockpile 

46 

329 

0.72 

1 

1 

47 

314 

0.65 

1 

1 

48 

312 

0.64 

1 

1 

49 

298 

0.58 

1 

1 

50 

284 

0.51 

1 

1 

51 

279 

0.48 

1 

1 

52 

276 

0.46 

1 

1 

53 

272 

0.44 

1 

1 

54 

270 

0.43 

1 

1 

55 

267 

0.42 

1 

1 

56 

262 

0.39 

1 

1 

57 

242 

0.28 

1 

1 

58 

237 

0.24 

1 

1 

59 

237 

0.24 

1 

1 

60 

232 

0.21 

1 

1 

61 

224 

0.16 

1 

1 

62 

222 

0.15 

1 

1 

63 

219 

0.13 

1 

1 

64 

219 

0.13 

1 

1 

65 

218 

0.12 

1 

1 

66 

215 

0.10 

1 

1 

67 

206 

0.04 

1 

1 

68 

204 

0.03 

1 

1 

69 

204 

0.03 

1 

1 

70 

204 

0.03 

1 

1 

71 

203 

0.02 

1 

1 

72 

200 

0.00 

0 

0 

446  414 


TOTALS 


349 


Note  that  the  attacker's  optimal  strategy  is  not  deter¬ 
mined  here.  For  a  given  number  k  of  attacking  missiles,  his 
problem  is 


maximize  V(d,a) 
subject  to 
T 

E  a(  i)  <_  k, 

i=l 

where 

T  /  a(i)  _ 

V( d, a)  =  E  v ( i )  (  1  -  n  (1  -  (0.8)d(1’J)) 
1=1  \  j=l 


with  d  as  in  Table  4.  This  is  a  non-linear  "knapsack"  type 
problem  whose  solution  could  be  found  by  dynamic  programming. 
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PROOF  OF  THE  LEMMA 


Lemma :  Let  a  =  (a^)  and  b  =  (b^)  be  a  pair  of  positive, 

monotonically  nonincreasing  vectors  of  size  t  such  that 


(a) 

(b) 


s 


I  a 

i=l 


i 


> 


s 

y  b. 

.  -i  i 

1=1 . 


t 

I  v 

i=l  1 


for  s  =  1,2,...,  t-1  • 


Then  for  q  «  (0,1), 


t 

n 

i=l 


< 


t 

n 

i=l 


Proof :  We  first  show  that  there  exists  a  txt  doubly  stochastic 

matrix  M  such  that 


b  =  Ma. 


Let 


M. 


l’a2 

ai-bi 


ai-bi 


Va2 


o  ...  o\ 
o  . . .  o 
l  . . .  o 


NO  0  0  ...  is 

stochastic.  Let  a"*-  =  M^a  and  note  that 

A-l 


Note  that  is  doubly 


>  b 


2 


In  general ,  we  will  have  a  vector  a^  b  such  that 


k-1  -  h 
ai  "  bi 


for  i  =  1,. . k-1 


3  k-1  3 

> 


l  =  I  >  I  b 

i-1  i=l  1  1=1  1 


for  s  =  kj . . . j t-1 


t  .  t 

i  4  -  i  \ 

i-i  1  1=1  1 


and  if  we  pre-multiply  such  a  vector  by  the  doubly  stochastic 
matrix 


we  obtain  a  new  vector  a  which  has  the  above  properties^ until 
k  =  t  -  1,  when  ab-b  =  b. 


A-2 


Since  the  product  of  doubly  stochastic  matrices  remains 

#  doubly  stochastic,  the  matrix 

M  =  Mt_1  . . .  Mx 

« 

has  this  property,  and  also 

# 

b  =  Ma  . 


Let  M  =  (X^j)  so  that 


b  .  =  7  X .  .  a  .  , 


where 


'  T  X,  .  =  1,  X.  .  >  0  . 
j=l  1J  1J  - 

z 

Since  the  function  £n(l-q  )  is  strictly  concave  as  a 
function  of  z,  we  have 


b .  t  /  a  .\ 

£n(l-q  x)  >  l  X±  An ^1-q  Jj 

X. 


so  that 


t  /  a  \  ij 
=  £n  n  I 1-q  J  J  , 
J  =  1  \  / 


a.  \^ij 


z  d  .  z  z  /  a .  \ 

l  An (1-q  >  l  An  II  (l-q  J  ) 

i  =  l  i=l  .1  =  1  V  / 

Li1-*3) 


t  t 
=  In  n 

1=1  j=l 


X.  . 

ij 


=  to  l  (l-q^') 

j=l\  / 


l  *ii 

1  =  1  1J 


A-3 


1,  so  that 


But  since  M  is  doubly  'Stochastic, 


I  A 


ij 


which  completes  the  proof. 

We  note  the  dependency  on  log-concavity  in  this  proof. 
This  property  arises  in  several  instances,  and  seems  to  be  an 
essential  feature  of  the  problem.  See  also  Appendix  B. 
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ANOTHER  ALGORITHM  FOR  THE  ALL-INTEGER  VERSION 
OF  THE  PRIM-READ  MODEL 


1 .  Introduction 

Falk  [3]  has  given  an  algorithm  for  finding  an  optimal 
all-integer  solution  of  the  Prim-Read  allocation  model.  We 
present  here  another  such  algorithm.  This  algorithm  requires 
time  proportional  to  log  p,  where  p  is  the  price;  the  algorithm 
of  [3]  requires  time  proportional  to  p.  Thus,  the  present 
algorithm  is  probably  superior  for  hand  computation  at  least; 
for  machine  computation,  both  are  so  fast  that  there  is  pre¬ 
sumably  no  great  advantage.  However,  our  main  purpose  in 
presenting  another  algorithm  is  not  computation  but  insight 
about  the  validity  of  the  original  (non-integral)  Prim-Read 
model.  This  model  has  many  very  attractive  properties  (see 
[4]  ) ,  but  of  course  the  non-integral  assignments  produced  are 
not  feasible.  The  integral  assignments  produced  by  our 
algorithm  can  be  readily  compared  with  the  semi-integral  model 
of  [1] .  The  latter  model  can  in  turn  be  (somewhat  less  readily) 
compared  with  the  original  Prim-Read  model.  Our  ultimate  goal 
(among  others)  is  to  provide  easily  computed , guaranteed  bounds 
for  the  optimal  integer  solution  in  terms  of  the  Prim-Read 
solution.  This  note  is  the  third  (after  [1,2])  of  a  series 
aimed  at  that  goal. 

2 .  Characterizing  an  Optimal  Solution 

Before  giving  the  algorithm  we  must  establish  some  results 
which  will  characterize  the  particular  solution  it  will  find. 

As  usual,  the  optimization  problem  to  be  solved  is  as  follows: 

(2.1)  minimize  D  =  J  d(i)  , 

i 

subject  to  the  constraints 


(2.2) 


V(d, a)  =  1  - 


a 

n  (l- 
i=l 


qd(1)) 


a/p 


for  all  positive  integers  a,  where  the  d(i)  are  required  to  be 

non-negative  integers.  As  discussed  in  [1]  and  [2].  p  need 

not  be  an  integer;  of  course,  p  >  0.  By  Proposition  3-1  of  [1],  # 

there  always  exists  a  solution  for  which  d ( J  )  >_  d(j+l)  for  all  j, 

that  is,  a  monotone  solution. 

We  may  describe  such  a  solution,  as  in  [1],  by  a  sequence 
of  positive  integers  x^ 


X2  ± 


—  xz 


(2.3) 


d(i )  =  2,  ,  0  <  i  <_  x^  5 
d(i)  =  &-1,  x£  <  i  £  x£-1  , 


d  ( i ) 
d(i) 
d  (i ) 


2  ,  x3  <  i  <  x2  , 

1  >  x2  <  i  £  , 

0  ,  x1  <  i  . 


Of  course,  it  is  clear  that  at  the  optimum,  x^  =  {p}  -  1,  % 

where  { p }  denotes  the  least  integer  >_  p  .  We  also  note  that 
D  is  given  by  x^  +...+  x. . 

Clearly,  we  could  now  restate  the  optimization  problem 
(2.1-2)  in  terms  of  the  x,  ,  and  indeed  this  was  done,  for  non-  ® 

integral  x^,  as  formulas  (4.1-4)  of  [1].  We  do  not  do  this 
now;  but  we  do  need  the  notion  of  an  x,  not  being  tightly  con¬ 
strained.  Call  x^  in  a  feasible  set  of  x^  free  if  the  set 

(x\  },  defined  by  x'.  =  x.  for  £  >  j  >  k,  x!  =  x.  -  1,  and  ® 

J  JJ  -  k  k 

xT.  =  x-,  for  k  >  j  >_  1 ,  is  also  feasible.  (Note  that  this 

J 

entails  x^  >  x^+^.)  That  is,  the  set  of  d's  starts  out  the 

same,  but  the  value  k  -  1  starts  one  step  earlier,  and  then  ® 


B-2 


continues  to  the  end.  If  x^  =  1,  the  effect  of  setting  = 
x  -  1  is  the  same  as  that  of  replacing  £  by  £  -  1. 

We  will  show  that  we  can  always  find  an  optimal  solution 
in  which  no  x^  is  free;  we  begin  with  a  simple  algebraic 
proposition.  We  note  in  passing  that  it  is  equivalent  to  the 
fact  that  1  -  is  log-concave  in  t . 

Proposition  2.1:  If  d  >  e  +  1  and  e  >_  0 ,  then 
(l-qd)(l-qe)  <  (l-qd_1) (l-qe+1)  . 


Proof : 


(l-qd  1 ) (l-qe+1 )  -  (l-qd)(l-qe) 


d-1  e+1  .  d  ,  e 

=  -  q  -q  +q+q 


e,..  d-e  d-e-1) 
=  q  (1-q+q  -q 


=  qe(l-q)(l-qd  e_1)  >  0 


0 


Armed  with  this,  we  can  prove  a  crucial  fact.  We  have  to 
state  the  result  a  bit  more  technically  than  one  would  like 
because  of  the  possibility  that  xk_^  =  x^_2  case,  no 

d(i)  is  equal  to  k  -  2). 


Proposition  2.2:  Let  x^,  ...,  x^  define  an  optimal 

(feasible)  solution  to  (2.1-2),  as  in  (2.3),  and  suppose  that 
x^  is  free,  where  k  >  1.  Let  k'  be  the  largest  value  for  which 
xk,  =  xk_2_-  (Usually,  k'  =  k  -  1).  Then  if  x^  is  replaced  by 
xk  -  1  and  xk,  is  replaced  by  x,  ,+1,  the  resulting  sequence 
(xl )  is  again  an  optimal  solution. 


B-3 


Proof :  Since  D  is  unchanged,  we  need  only  check  that  the 


resulting  set  of  x.  is  feasible.  Consider  first  the  (non- 

monotone)  sequence  of  d(i)  that  comes  from  (2.3)  by  setting 

d(xk~l)  =  k  -  1  and  d(xk_^)  =  Clearly  D  is  unchanged. 

Moreover,  this  new  solution  obviously  satisfies  (2.2)  for 

a  _>  and  it  satisfies  (2.2)  for  a  <  xk-1  because  x^  was 

free.  Now  we  make  this  allocation  monotone.  Set  d(x,  ,)  = 

k-1 

k  -  1,  d(xk_1+l)  =  k'.  By  the  definition  of  k*  ,  xk_^  +  1  = 
x^,  +  1  and  k'  is  one  greater  than  the  old  value  of 
d(xk_-^+l)  =  d(xk,  +1 ) ,  so  that  the  new  assignment  has  the  same 
stockpile  D  as  before.  But  this  assignment  is  also  feasible 
by  Proposition  2.1.  Since  this  is  just  the  assignment  speci¬ 
fied  by  the  sequence  {x! },  the  proof  is  complete.  0 

The  above  proposition  leads  immediately  to  the  following 
theorem,  which  completely  characterizes  one  solution  to  the 
all-integer  problem. 


Theorem  2.1:  There  is  a  unique  feasible  allocation  of  typ 
(2.3)  for  which  no  x^  is  free,  and  this  allocation  is  optimal. 

Proof :  Clearly,  there  is  only  one  choice  of  x^  that  is 

(feasible,  but)  not  free.  Once  x^  is  specified,  there  is  only 

one  choice  of  x^_^  that  is  not  free,  and  so  on.  Therefore, 

the  x.  are  completely  specified  in  turn  by  the  conditions  on 

them.  (We  will  shortly  consider  how  to  calculate  the  x..)  To 

<3 

see  that  this  solution  is  optimal,  we  need  only  show  that  an 
optimal  solution  can  be  converted  to  it  while  remaining  optimal 


Consider  any  such  optimal  solution.  If  no  x.  is  free, 
we  are  done.  Let  k  be  the  largest  index  for  which  x^  is  free. 
We  cannot  have  k  =  1,  for  then  the  solution  would  not  have 
been  optimal.  By  Proposition  2.2,  it  is  possible  to  reduce 
x^  by  one  and  still  have  an  optimal  solution.  If  the  resulting 
x^  is  still  free,  we  may  continue  the  process  until  it  is  no 


longer  free.  (If  k  =  l,  this  could  involve  reducing  x^  to  zero 

and  hence  replacing  Z  by  Z  -  1.  )  In  any  case,  this  process 

leads  to  an  optimal  solution  in  which  the  largest  j  for  which 

x j  is  free  is  <  k.  If  we  continue  to  do  this,  we  eventually 

arrive  at  an  optimal  assignment  for  which  no  x.  is  free.  This 

1 

completes  the  proof.  Q 


3 .  The  A1 qori thm 

Theorem  2.1  characterizes  an  optimal  solution,  and  one 
may  clearly  derive  an  algorithm  from  it.  Like  Falk's  algorithm 
[3],  it  is  a  greedy  algorithm;  but  it  is  greedy  in  the  x.  of 
(2.3),  not  in  the  d(i)  of  (2.1-2).  Its  solution  will  in 
general  be  different  from  Falk's,  since  it  produces  monotone 
d(i),  but  of  course  the  resulting  stockpile  D  will  always  be 
the  same. 

Of  course.  Theorem  2.1  does  not  give  an  algorithm  in  itself, 
since  we  must  make  explicit  how  to  choose  the  x^  so  that  they 
are  not  free.  We  now  state  a  complete  algorithm.  Although 
it  looks  tedious,  the  algorithm  is  actually  rather  easy  to 
apply.  We  use  the  notations  [t]  and  {t}  to  denote  the  greatest 
integer  <_  t  and  the  least  integer  >_  t,  respectively. 


Algorithm:  Set  l  =  {-log  p/log  q}  and  set  6.  =  -log(l-qJ), 

1  =1,  2,  ...,  1.  Also,  set  x?+1  =0,  0£  =  0.  Then,  for  j  =  l, 


Z  -  1, 

(3.1) 

(3.2) 

(3-3) 


,  2  set  in  turn 


a  . 

i-i 

♦i-i  -  «* 


h'-i 


-  17  j-i  > 

j-iCaj-i]  +  log(1-[aj-i]/p)  • 

g j _t_ { a ^  log(l—  p )  if  <  1  , 

0  otherwise. 


B-5 


(3-4) 


(3-5) 


(  max(<f> ’  ,  cj)  ”  _  )  -  9  . 

x .  =  { - 4 - ^ — j - - 

J  5,  n  -  5, 


j-1 


9.  -i  =  0.  +  (6.  , -6 . )x. 
0-1  J  J-1  J  J 


Finally,  set  x^  =  { p }  -  1.  Note  that  the  calculation  of  9n  may 
be  omitted. 

We  must  show  that  this  works;  we  state  that  as  a  theorem. 


Theorem  3.1:  The  above  algorithm  produces  a  feasible 
;nce  of  x.  such  that  noi 
sent  an  optimal  allocation. 


sequence  of  x.  such  that  none  are  free;  they  therefore  repre- 


Proof :  Of  course,  £  must  be  great  enough  so  that  the 

solution  is  feasible  for  the  first  attacker,  so  that  we  must 
have ,  from  (2.2), 

,d(  1 ) 


1/p  1  1  -  ( 1-q 


=  q 


so  that 


-  log  P  2l  &  l°g  q5 


and  hence 

£  _>  (-log  p/log  q)  . 

On  the  other  hand,  if  £  is  strictly  greater  than  this,  x^  must 
be  free  since  even  x^  =  0 ,  =  { P }  -  1  is  feasible.  Hence, 

£  =  (-log  p/log  q} . 

The  remainder  of  our  proof  is  closely  related  to  ideas 

developed  in  Section  4  of  [1]  .  In  particular,  V(.a,d)  for  any  a 

in  the  range  x.  <  a  <_  x._,  is  given,  for  some  0,  by 
J  J 


B-6 


since 


V  ( a ,  d )  =  1  -  exp(0-6._,a)  , 

J  -L 


X  . 

J 


a-x . 


V(a,d)  =  1  -  HJ(l-qd(l))  •  (l-qJ_1)"  ^ 

i=l 


X  . 


1  -  H  (l-qd^b  •  (1-q^  d)  ^*exp(-6.  ,a) 
1=1  J-1 


where  we  take  0  to  be  the  natural  logarithm  of  the  product  and 
the  following  factor.  We  will  set  0  =  0._,  in  any  such  range; 
we  must  show  that  0._,  is  as  defined  in  the  algorithm.  Of 
course ,  0^  =  0  in  the  first  range,  as  it  should  be.  For  any 
other  range,  x.  is  the  value  of  a  at  the  intersection  of  the 

J 

curves  v  =  1  -  exp ( 0 . -6 . a)  and  v  =  1  -  exp(9.-6.a);  this 

J  -1-  J  -*•  J  J 

easily  leads  to 


(3-6) 


a  =  x .  = 

J 


0  .  ,-0 . 

J"1  J 
6  .  ,-<S.  5 
0-1  J 


which  is  equivalent  to  (3-5).  (Formula  (3.6)  already  occurs 

in  the  proof  of  Theorem  4.1  of  [1],  except  that  there  it  was 

stated  in  terms  of  y,  =  -  0..) 

J  J 

Of  course,  before  we  can  determine  0.  we  must  determine 

J  -L 

x. .  We'  must  choose  x.  as  small  as  possible  consistent  with  the 
J  J 

requirement  that  the  curve  v  =  1  -  exp  ( 0  .  - -6  .  .  a )  lies  on  or 

J  -i  J  -L 

below  the  line  v  =  a/p  at  all  integer  a.  If  they  were  tangent, 
the  point  of  tangency  would  be  at  a.  =  p-1/6.  _  •  this  is 
formula  (4.5)  of  [1],  However,  since  only  integer  values  of 

a  are  to  be  considered,  we  need  that  the  curve  lies  on  or  below 
the  line  at  a  =  [aj_^]  and  a  =  (a.  p  ,  unless  (a._p  >  p,  in 

which  case  this  point  should  be  ignored.  If  the  curve  v  =  1  - 


exp(0-6  a)  goes  through  the  point  (a,v)  =  (a, a/p),  then 


we 


have 


exp (0-6  a)  =  1  -  a/p 
J 
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so  that 


(3-7)  0  =  6.  -.a  +  log(l-a/p)  . 

J 

From  this  we  see  that  <f>'  and  <j>"  in  (3.2-3)  are  the  0 

J  J  a 

of  (3.7)  corresponding  to  a  =  [a.-,]  and  a  =  {a.,},  respectively, 

provided  (a,  ,}  <  p.  Consequently  0  =  max(<j>!  ,(p"  ,  )  is  the 
J  ^  J  ~  -L  J  —  1 

smallest  0  for  which  v  =  1  -  exp(0-S.  .a)  does  not  go  above 

v  =  a/p  at  any  integer  point.  (Note  that,  in  the  case  {a.  , } 

J  J- 

>,  P5  9  =  4>j_2  since  the  latter  is  positive,  which  is  as  it 
should  be.)  Hence,  by  (3.6),  we  need  to  have,  and  need  only  have, 

max(<j>'_  1  ,<fr"  )  -  0. 

Y  .  _ .  J  J  j _ 


in  order  for  x.  to  produce  a  feasible  solution.  Therefore,  if 
J 

we  take  x.  to  be  as  small  as  possible,  but  still  an  integer, 
x.  will  not  be  free.  This  gives  us  (3.4).  Finally,  it  is 

J 

clear  that  x^  =  { p }  -  1.  This  completes  the  proof.  [[ 

The  form  of  (3-1-5)  permits  rather  ready  comparison  with 
the  results  of  the  algorithm  for  the  semi-integral  model 
presented  in  Theorem  4.2  of  [1] .  This  fact  is  a  primary 
motivation  for  developing  the  above  algorithm;  however,  such 
a  comparison  will  be  deferred  to  a  later  note. 


4 .  Numerical  Examples 

We  discuss  here  two  numerical  examples,  hand-calculated. 
In  a  later  note  we  will  give  more  extensive  numerical  results 
produced  with  the  aid  of  a  computer.  First,  consider  the 
case  q  =  .25,  P  =  50,  considered  in  [1]  and  [3] •  We  have 
i  =  3,  and  ^  =  .2877,  «2  =  .06454,  §3  =  .01575,  x4  =  0, 

0^  =  0.  Then, 

a2  =  34.51  , 

$'2  =  1-0549,  =  1.0549  , 
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{21.62}  =  22 
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02  =  1.0734  , 
ax  =  46.52  , 

*j_  =  10.708,  <j>!J  =  10.708  , 
x2  =  -{43.17}  =  44  , 

and,  finally,  x1  =  49.  The  total  stockpile  is  D  =  22  +  44  +  49  = 
115. 

Some  remarks  are  in  order  here.  First,  we  have  left  off  the 
unnecessary  calculation  of  Second,  the  stockpile  agrees 

with  the  corresponding  one  in  [3],  as  it  must.  Moreover,  the 
numbers  of  3's,  2's,  and  l's  tally  exactly  with  those  in  [3], 
but  of  course  -the  sequence  is  different.  It  seems  likely  that 
such  an  exact  tally  should  always  occur,  but  this  has  not  been 
proved.  Finally,  note  the  surprising  fact  that  in  each  case, 

<)>'.  =  <j>'.’ .  This  is  not  a  coincidence;  in  fact  we  will  prove  in 
the  last  section  that  it  always  happens  when  p  -  q-^  is  an 
integer.  Thus,  when  q  =  .5,  .25,  .2,  .1,  ...  and  p  is  an 
integer,  this  applies,  a  useful  fact  for  hand  calculation. 

Now  we  let  q  =  .2,  p  =  50,  again  considered  in  [1]  and 
[3].  We  have  £  =  3,  S-L  =  .2231,  62  =  .04082,  5,  =  .008032, 
x2|  =  0,  =  0.  Then, 

a. 2  =  25.50  , 

*'2  =  <l>2  "  -3274  , 

x3  =  {9.985}  =  10  , 

e2  =  .3279  , 

a1  =  45.51  , 

=  $1  =  7.737  , 
x2  =  {40.65}  =  4l  , 

and  x1  =  49.  This  time,  D  =  10  +  41  +  49  =  100. 
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As  before,  even  the  individual  Xj  agree  exactly  with 
the  corresponding  numbers  of  3's,  2's,  and  l's  in  [3], 


5 .  Conditions  for  a  Simplified  Calculation 

It  has  been  already  said  that  <p  '.  =  <j/.'  whenever  p  -  q 

J  J 

is  an  integer.  Although  this  is  a  rare  event  for  p  and  q 
completely  arbitrary,  it  will  often  occur  for  the  values 
typically  chosen,  simplifying  hand  calculation.  'For  this 
reason,  it  is  worthwhile  to  prove  this,  and  a  little  more, 
explicitly . 


-j 


Proposition  5.1:  For  1  <  j  <  I,  it  will  happen  that 

4> I  =  in  the  algorithm  of  Section  3  if  and  only  if  either 

P  -  1/6 .  or  p  -  q  d  is  an  integer. 

J 


Proof :  If  p  -  1/6.  =  a.  is  an  integer,  then  [a.] 

J  J  J 

and  the  result  is  trivial.  Otherwise,  we  may  set  [a^.] 

{a.}  =  a  +  1.  In  this  case,  <f> !  =  <p  yields 
J  J  J 


{a  . } 

J 


0  -  <p -  <p '.  -  6.(a+l)  +  log(l  -  (a+l)/p)  -  6.  a  log  (1-a/p) 
J  J  J  J 

=  5,  +  log  H-.~.  ,a  ~  1  _  i0g  ~  a 
J  P  p 

„  ,  n  p-a-1 
j  &  p  -  a  5 

so  that 


t.  n  p-a-1 
j  &  p  -  a 


But  6.  =  -log(l-q^),  so 
J 


i  -  qJ  -  ^---a  -  1  =  1  — , 

p  -  a  p  -  a  5 
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which  is  equivalent  to 


a  =  p  -  q  . 

Since  a  is  an  integer  it  is  necessary  that  p  -  q-^  be  an 
integer  in  this  case. 

To  prove  that  this  condition  is  also  sufficient,  we 
need  only  show  that 

p-q^^_p-l/<5.<_p-q^  +1, 

J 

which  is  equivalent  to 

qJ’  1  <  l/(-q_J+l)  =  qJ/(l-qJ')  . 

But  since  0  <  q^  <  1,  we  may  expand  6.  =  -log(l-q^)  and 
*?  -?  J 

qJ/(l-qJ)  in  series;  the  above  is  then  equivalent  to 

qJ’  <  qJ’  +  q2J/2  +  q3j'/3  +  ... 

<_  qJ'  +  q2j"  +  q3j'  +  .  .  .  ; 

this  is  obvious,  so  the  proof  is  complete.  f] 
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